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1 Introduction 

Rigid body dynamics in an ideal incompressible fluid is rich with problems interesting from 
mathematical point of view, in particular, the research of integrable problems. Certainly, the 
most famous are the three integrable cases under the names of Kirchhoff, Clebsch, Steklov and 
Lyapunov. The latter two cases are more interesting because there are no a obvious symmetry 
groups associated with the additional integrals of motion. 

All these classical cases were discovered and carefully studied in the 18th and 19th cen- 
turies [5]. For instance, the Kirchhoff equations for the Clebsch and Steklov-Lyapunov cases 
were first solved explicitly by Kotter after some mysterious separation of variables [11, 12]. At 
the moment no separation which is alternative to his original separation of variables is known 
for these systems, even though there is a lot of literature dedicated to the problem, including 
the theta-functions solutions, associated with the Lax matrices, and the detailed geometric de- 
scriptions of the invariant surfaces on which the motions evolve, see books [1, 2] and references 
within. 

In this paper we apply the bi-Hamiltonian geometry to the direct calculation of variables of 
new separation for the Steklov-Lyapunov system and it's Rubanovsky generalization. We have 
to notice right away, that our main purpose is the development of the bi-Hamiltonian geometry 
instead of integration of particular equations of motion using separation of variables method. 

An integrable system is separable, if there are n separation relations 



^i{ui,pui,Hi,...,Hn) = 0, i = l,...,n, with det 



/O, (1.1) 



connecting single pairs {ui,puj of canonical variables of separation with the n functionally 
independent Hamiltonians Hi , . . . , Hn ■ Solving these relations in terms of . one gets the Jacobi 
equations and the corresponding additively separable complete integral of the Hamilton-Jacobi 
equation 



= ^ / Pu,{u'i, ai,..., an)du'i, aj = Hj. 
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Of course, any variables of separation are determined up to the trivial transformation 

Ui ^ Ui = fi{ui,pui), (1.2) 

which preserve all the properties of algebraic curves defined by separated relations (1.1). Howe- 
ver, if we have another separated relations for the same integrals of motion 

^i{vi,pvi,Hi, . . . ,Hn) = 0, i = l,...,n, with det 

which can not be reduced to initial ones (1.1) by trivial change of variables (1.2), we usually 
say about different variables of separation. Of course, any two families of canonical variables on 
a given phase space are related by a generic canonical transformations 

Vi = gi{ui, . . . ,Un,Pl, ■ ■ ■ ,Pn), Pvi = hi{ui, . . . ,Un,Pl, ■ ■ ■ ,Pn)- (1.3) 

In contrast with (1.2) the notion of generic transformations (1.3) allows us to study relations 
between distinct algebraic curves, for instance covering of the algebraic curves (see works of 
Poincare, Humbert, Frey, Kani, Kuhn, Shaska) or curves with isogenous Abel varietes (see 
works of Richelot, Brock, Hayashida, Nishi, Ibukiyama, Katsura, van Wamelen). Such relations 
give us a lot of examples of the reductions of Abelian integrals (see works of Hermite, Goursat, 
Burkhardt, Brioschi, Bolza) and, therefore, they may be a source of new ideas in number theo- 
ry, algebraic geometry and modern cryptography (see works of Tate, Faltings, Zarhin, Lange, 
McMullen, Merel). 

Nevertheless, motivation for the search of such different variables of separation do not pure 
mathematical, because different variables of separation may be useful in different perturbation 
theories [18, 22, 24], as well as in distinct procedures of quantization and various methods of 
qualitative analysis, etc. 

The milestones of the variables separation technique include the works of Stackel, Levi- 
Civita, Eisenhart, Benenti and others. The majority of results was obtained for a very special 
class of integrable systems, important from the physical point of view, namely for the natural 
Hamiltonian systems with quadratic in momenta integrals of motion on cotangent bundles to 
Riemannian manifolds. The Kowalevski and Chaplygin results on separation of variables for 
the systems with higher-order integrals of motion have been missed out of this scheme until 
recently [18, 22, 24]. 

In the Steklov-Lyapunov case we have quadratic integrals of motions, but the phase space 
is the Poisson manifold instead of the cotangent bundle to Riemannian manifold. So, in this 
case we can use neither the Levi-Civita criteria, nor the Eisenhart-Benenti theory. Below we 
show how variables of separation for the given integrable system may be calculated without any 
additional information (Killing tensors. Lax matrices, r-matrices, links with soliton equations 
etc.). 



2 Steklov— Lyapunov system 

Following Kirchhoff, we consider the potential motion of a finite rigid body submerged in an 
infinitely large volume of irrotational, incompressible, inviscid fluid that is at rest at infinity, so 
that the induced motion of particles of the fluid is completely determined by the motion of the 
body [9]. In this case, the motion of rigid body is described by the classical Kirchhoff equations 

M = Mxil + pxC/, p = pxQ, (2.1) 

here x x y stands for the vector product of three-dimensional vectors. Vectors M and p are 
the impulsive momentum and the impulsive force while O and U are the angular and linear 
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velocities of the body. All these vectors in are expressed in the body frame attached to the 
body originating at the center of buoyancy [9] . 

A rigid body in the ideal fluid is an important example of Hamiltonian systems on a dual to 
Lie algebra e(3) = so(3) >< M^. The dual space e*(3) is the Poisson manifold endowed with the 
canonical Lie-Poisson brackets 



{Mi,Mj} = EijkMk, {Mi,pj] = EijkPk, {Pi,Pj} = 0, 
where Sijk is a totally skew-symmetric tensor. There are two Casimir elements 

3 3 

Ci = {p,p) = \p\^ = Y.pI ^2 = {p,M) = ^p.M„ 



i=l 



1=1 



(2.2) 



(2.3) 



usmg 



(2.4) 



where x is a cross product, [•, •] is a matrix commutator and index n means a 3 x 3 antisymmetric 
matrix associated with the vector z. Using this agreement we can rewrite canonical Poisson 
bivector on e* (3) in the following compact form 



where (x^y) means scalar product of two three-dimensional vectors x,y G M^. 

As usual [9], element M E so(3) is identified with three-dimensional vector M G 
well known isomorphism of the Lie algebras (M^, x) and so(3), [•, •] 
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(2.5) 



The Hamilton function H = H{p,M) and Lie-Poisson brackets (2.2) allow us to define the 
Hamiltonian equations of motion 



M = Mx— +PX-, 



p = p X 



dH 
dM' 



(2.6) 



These generic Euler's equations on e*(3) coincide with the Kirchhoff equations (2.1), if H{p, M) 
is a second-order polynomial in variables M and p. 

Remark 1. The Lie-Poisson dynamics on e*(3) can be interpreted as resulting from reduction by 
the symmetry group -^(3) of the full dynamics on the twelve-dimensional phase space T*E{3). 
Here, symmetry means that the Hamiltonian that describes the dynamics in T*E{3) is an 
invariant to actions of -£^(3), i.e., one can translate the inertial frame or rotate it in any direction 
without affecting the equations of motion [15]. 

The Steklov-Lyapunov case of the rigid body motion is characterized by the following second- 
order homogeneous polynomial integrals of motion 



Hi = {M,M) - 2{Ap,M) - ((A2 + 2A^)p,p), 

H2 = {AM,M) + 2{A^p,M) - ((A^p,^) -trA2(Ap,p)), 



(2.7) 



4 



A.V. Tsiganov 



where wedge denotes an adjoint matrix, i.e. a cofactor matrix = (det A)A ^. For integra- 
bility A has to be symmetric matrix, which may be reduced to the diagonal form 

A = diag(ai, 02, as), Oj G M, (2.8) 

using hnear canonical transformations of e*(3). From physical point of view it means that the 
body axes can always be chosen so that A is diagonal. 

Remark 2. In [17] Steklov found integrable Hamiltonian H2, whereas Lyapunov proved inte- 
grability of the Kirchhoff equations with Hamiltonian Hi in [14]. The generic family of the 
Steklov-Lyapunov integrable systems was studied by Kolosov in [10]. 




2.1 Separation of variables by Kotter 

The explicit integration of the classical Steklov-Lyapunov systems via separation of variables 
had been first made by F. Kotter in 1900 [11]. Here we want to add some new details in the 
known coincidence of the Kotter variables of separation vi^2 with the elliptic coordinates on the 
sphere. 

According to [23, 25] there is a Poisson map, which identifies the Steklov-Lyapunov system 
with a system that describes motion on the surface of a unit two-dimensional sphere in 
a fourth-degree polynomial potential field. This dynamical system is separable in standard 
elliptic coordinates on the sphere, and the inverse Poisson map allows us to get complete solution 
of the Steklov-Lyapunov system. 

Proposition 1. 7/B = tr A — A and C = — 4A^ are diagonal matrices with entries 

'a2 + ag 
B = I ai + 03 

ai + a^f 

then the Poisson map f : (p, M) — )• (x, J), defined by 

_ (M - Bp) X p 
^ ~ \{M - Bp) X p\ ' 

where 

n=3 

j,k=i 

relates manifold e*(3) with coordinates {p,M) and cotangent bundle T*S^ to the unit two- 
dimensional sphere with coordinates {x, J). 

The proof consists in the verification of the Lie-Poisson brackets between variables x and J 

and calculation of the corresponding Casimir functions 
(x, x) = \x\^ = 1, and (x, J) = 0. 



J = M + C[x,xxp] (2.9) 



So, this coajoint orbit of e*(3) with coordinates x and J is simplectomorphic to cotangent 
bundle T*S^ to the unit two-dimensional sphere S^, see for instance [15]. 
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Inverse Poisson map f-^ : (x, J) {p, M) looks like 

p = aJ + (5{xx J), M = J -C[x,xxp\^, (2.10) 
where functions q, /3 on x, J are solutions of the following equations 
Ci = = aVl^ + X 

C2 = {p,M) = a'^{\J\'^{x,Ax) - {{x X J), A(x x J))) +2a/3(J, A(x x J)) 
+ /32(((x X J),B(x X J)) - 2(J, AJ)) +a|J|2. 

Proposition 2. The Poisson map (2.9) relates the Steklov integral of motion Hi{p,M) (2.7) 
with the natural Hamilton function on T*§^ 

Hi{x,J) = (J, J) +4(x,Bx)((trACi - C2) - Bx)) +4Ci(x, A^x) 

+ 2trAC2 -tr(A2 + 4A^)C7i. (2.11) 

The Lyapunov integral H2{p,M) (2.7) is equal to 

H2{x,J) = (J, A J) -4(x, A^x)(C7i(x,Bx) - (trACi -C2)) 

+ 2trAC2 -2detACi. (2.12) 

Here Ci^2 o-f^ values of the Casimir functions (2.3). 

The proof of this proposition and all the details may be found in [23, 25]. 

Remark 3. It is well-known that the generic level sets of the Casimir functions (coadjoint orbits) 
on e*(3) are only diffeomorphic to the cotangent bundle r*S^. It allows us to directly connect 
the Kirchhoff equations with the equations of motion by geodesic on S^, however, at C2 7^ 
we have to destroy the standard symplectic structure on T*S^ by adding some "monopole" 
terms [15]. 

In the Steklov-Lyapunov case we use the Poisson map (2.9), which preserves the standard 
symplectic structure on r*S^. As a punishment for this preservation of the standard symplectic 
structure we have to consider potential motion on instead of the geodesic motion. 

The Hamiltonians Hi^2ix, J) (2.11), (2.12) on T*S^ are separable in the elliptic (spheroconi- 
cal) coordinates ^1^21 which are zeroes of the function 

This variables of separation satisfy the following separated relations 

^{vi,pyJ = det{v,-A)pl + {£{v,),£{vi)) = 0, i = l,2, (2.14) 
where the three-dimensional vector £{X) is the so-called Kotter vector with entries 



^i(A) = ^^^^ {Mi + (2A + ai - tr A)pi) 

depending on the auxiliary variable A (spectral parameter). The explicit solution of the cor- 
responding Abel-Jacobi equations in theta- functions was given in [11]. In order get initial 
variables p, M as functions on time variable we have to substitute solutions of the Abel-Jacobi 
equations into the variables on T*S^ 



Uvi - ai){v2 - Oi) 2eijkXjXk{aj - ak) (, . , . x 

Xi = W 7 TT r, Ji = ((flj - Vi)p^^ - [Oi - V2)Pv2), 

\] [aj - ai){ak - Oi) vi - V2 ^ 

where {i,j,k) is permutation of (1,2,3), and then into the variables p and M (2.10). Another 
modern verification of the Kotters calculations may be found in [4, 8]. 
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3 Calculation of the variables of separation 
in bi-Hamiltonian geometry 

We can only guess how Kotter invented the variables of separation vi^2, which coincide with the 
elliptic coordinates on the auxiliary two-dimensional sphere, because he gave no explanations 
of calculations in very brief communication [11]. It is clear that behind the striking formulas 
there must be a certain geometric idea, but the domain of applicability of this idea is usually 
restricted by a partial model under consideration. For instance, we can not apply the Kotter 
separation to the Kowalevski top and Kowalevski separation to the Steklov-Lyapunov system 
etc. 

Our aim is to discuss some algorithm of calculation of the variables of separation in the 
framework of the bi-Hamiltonian geometry, which is applicable to many known integrable sys- 
tems [7, 18, 20, 21, 22, 24, 26]. In fact this algorithm consists of the following steps: 

• calculate the second Poisson bracket compatible with canonical one starting with the given 
integrals of motion in the involution with respect to this canonical Poisson bracket; 

• if the Poisson brackets have different symplectic leaves, calculate a projection of the second 
bracket on symplectic leaves of the first bracket; 

• calculate coordinates of separation as eigenvalues of the corresponding recursion operator; 

• calculate the canonically conjugated momenta with respect to the first Poisson bracket; 

• calculate the separated relations. 

The input of algorithm is a set of integrals of motion and canonical Poisson bracket, whereas 
output is a set of separated relations. All details about construction of a suitable projection are 
discussed in [7]. 

Because variables of separation (1.1) are defined up to canonical transformations Ui — )• 
f{ui,Pui) on the first step of this algorithm we have to narrow the search space using some 
artificial tricks. It is a main technical problem of this method. The second technical problem is 
the calculation of the momenta conjugated to obtained coordinates, see [7, 18, 22, 24]. 

3.1 Polynomial and rational Poisson brackets on e*(3) 

Bi-Hamiltonian structures can be seen as a dual formulation of integrability and separability, 
in the sense that they substitute a hierarchy of compatible Poisson structures to the hierarchy 
of functions in involution, which may be treated either as integrals of motion or as variables of 
separation. So, our first step is calculation of the second Poisson bivector P' compatible with 
kinematic Poisson bivector P. 

According to [20, 26] any separable system is a bi-integrable system, i.e. integrals of mo- 
tion (1-1) are in bi- involution 

{Hi,Hk} = {Hi,Hky = 0v, i,k = l,...,n, (3.1) 

with respect to compatible Poisson brackets {•, •} and {•,•}' associated with the Poisson bivec- 
tors P and P' , so that 

[P,P1 = 0, [P,P'1=0, [P',P'1 = 0. (3.2) 

Here [•, •] is the Schouten bracket. The definition of the second bracket {•,•}' in term of variables 
of separation may be found in [20, 26]. 

For the given integrable system fixed by a kinematic bivector P and a tuple of integrals of 
motion Hi , . . . , Hn bi-Hamiltonian construction of variables of separation consists in a direct 
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solution of the equations (3.1) and (3.2) with respect to an unknown bivector P' . The main 
problem is that the geometrically invariant equations (3.1), (3.2) have a'priory infinite number 
of solutions [18, 20, 21, 22, 24, 26]. 

In order to get a search algorithm of effectively computable solutions we have to narrow the 
search space by using some non- invariant additional assumptions. According to [20, 21] hereafter 
we assume that P' has the same foliations by symplectic leaves as P, i.e. that 

p'dCi,2 = (3.3) 

and P' doesn't have any other Casimir elements. The geometric meaning of this restriction is 
discussed in [20, 21]. In fact it allows us to avoid calculations of the projection of the second 
bracket on the symplectic leaves of the first bracket. 

In the Steklov-Lyapunov case solving equations (3.1), (3.2) and (3.3) in the space of ho- 
mogeneous second-order polynomial bivectors and of rational bivectors with the second-order 
homogeneous numerators and linear denominators we obtain the following two propositions. 

Proposition 3. // c and d are two numeric three-dimensional vectors, so that (c, c) = 0, then 
equations (3.2) and (3.3) on e*(3) have a polynomial solution 



Pi 



I {c,p)p^ (c, M)p^ + ipx M)0c+l{^+ a{c, d)) {p^p-{p,p))\ 

+ a(c X p) {d X p) 

* (c, M)M^ + {d,p)p^ + + d)) {p X M)^ 

y - a((c xp)x{dx M))^ j 



(3.4) 



a G C, and a rational solution 

compatible to each other, i.e. |-Pi,i'2l ~ 0- 

As above x is a cross product, the antisymmetric matrix is defined by vector z (2.4) and 
the matrix [x ® y)ij = Xiyj is determined by a pair of vectors x and y. In order to explain this 
notations we write out the corresponding Poisson brackets 

{Pi,Pj}i = £ijk{c,p)Pk, 

{pi, Mj}[ = Eijkic, M)pk + {px M)iCj + ^ (1 + a^(c, d)) [piPj - ^ 



p1 



\ 1=1 / 

+ a{cx p)i{dx p)j, (3.6) 
{Mi, Mj]'^ = Eijkiic, M)Mk + {d,p)pk) 

+ eijk (1 + a^(c, d)) {p X M)k - a((c x p) x {d x M))^ . 

The second brackets are equal to 
{PhPj}'2 = £ijkPk, 

M,,'. ^^^^ ((. M) . ^^^^(^ X ^)) . (3.V) 
We have to stress that this brackets are defined over complex field because (c, c) = cf +C2+C3 = 0. 



A.V. Tsiganov 



Substituting the Poisson brackets {•,-}'i2 (3-6), (3.7) into (3.1) and solving the resulting 
equations in the space of the second-order homogeneous polynomials Hi^2 one gets the following 
proposition. 

Proposition 4. The Steklov-Lyapunov integrals of motion Hi^2ip,M) (2.7) satisfy the equa- 
tion (3.1) at a = 1 and 

c = ^, X (\/a2 - as, Vas - ai, Vai - ^2) , 

i/(ai - a2)[a2 - a3)[as - ai) 

d = (ax - a2){a2 - as){a3 - ai) ( . :, , , ^ ) 

Vv«2-a3 V«3-ai V«i-a2/ 

The proof is a straightforward calculation. 

Remark 4. In fact, polynomial bivector P[ has been obtained in [21] as an incidental re- 
sult by investigation of the Poisson bivectors on the Lie algebra so* (4) and the corresponding 
integrable cases in the Euler equations on so* (4). Now we recover this bivector by solving 
equations (3.1), (3.2) and (3.3) for the Steklov-Lyapunov system. 

Let us briefly discuss the bi-Hamiltonian structure related with the Kotter variables of sepa- 
ration. Elliptic coordinates on the sphere (2.13) are associated with the polynomial Poisson 
bivector 

Pin = ^xP, 

where Cx is a Lie derivative along the vector field X = Yl-^'^'^j with the following entries: 

X' = 0, X^+^ = [x X A{x X J)]., i = 1,2,3. 

Bivector is compatible with P and has the same foliation by symplectic leaves as P. 

Using the Poisson map (2.9), (2.10) we can easily express P'^i^ in the initial variables p, M. 
It will be a rational bivector P^^^ = R/Q, where i? is a bivector with fourth-order homoge- 
neous polynomial entries and Q = \{M — Bp) x is a fourth-order polynomial as well. So, 
we could directly calculate the Kotter variables solving equations (3.1), (3.2) and (3.3) in the 
corresponding space of rational bivectors. 



(3.8) 



3.2 Calculation of variables of separation 

The bi-involutivity of the integrals of motion (3.1) is equivalent to the existence of control 
matrix F defined by 

n 

P'dH = P{FdH) , or P'dHi = P ^ FijdHj, i = l,...,n. 

i=i 

The additional assumption (3.3) ensures that P is a non-degenerate matrix and the eigenvalues 
of F are the desired variables of separation [20, 21]. Moreover, for the so-called Stackel sepa- 
rable systems the suitable normalized left eigenvectors of the control matrix F form the Stackel 
matrix S [20, 21, 22, 24]. In this case separated relations (1.1) are affine equations in integrals 
of motion Hf.. 

Let us calculate the control matrices for the Steklov-Lyapunov system and, for brevity, 
introduce three constants 

3 

Tfc = trA'' EE ^of, A; = 0,1,2, 

i=l 
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and some linear functions on variables p, M 

Pk = {c,A^p), ak = {c,A^M), k = 0,1,2, (3.9) 
which are related to each other via the Casimir functions (2.3) on e*(3). For instance, 

3 1 

i=l 

In this notations control matrices associated with the Hamiltonians H12 (2.7) and the Poisson 
bivectors P[ 2 (3.4), (3.5) look like 




2/00 

-<7o - pi + npo. 



Po 3 
cri+ P2- npi 2ri 



(3.11) 



V PO 3 / 

Now we can simply calculate the desired variables of separation using two control matrices. 
Namely, let ui^2 be eigenvalues of the control matrix F2 

B{\) = det(F2 - A) = (A - ui)(A - ^2) 

\ Po 3 J 9 3po Po ' 

whereas the eigenvalues of Fi be doubled momenta 2pui 2 > so that the characteristic polynomial 
has the form 

A{X) = det(Fi - A) = (A - 2p,J(A - 2p„J 

= A^ + (ctq + 3pi - Tipo)A + 2/>i(cro + pi) - 2po{ai + P2). (3.13) 

Another equivalent definition of momenta p^- is given by a relation 

{c,p) {c,M) {c,dxp) . 
Pu, = ^^Ui 2 — ^ 3 ' ^ = 1'2. 

Now we can prove the following 

Proposition 5. On symplectic leaves o/e*(3) variables ui^2 o-ndpu-i2 ^'''^ canonical variables 

{Ui,Pui} = '^-, {'Ui-,Pui}'l=2pui, {Ui,p„.}2 = Uj, i = l,2. 

with respect to canonical Poisson bracket (2.2). 

The proof consists of the calculation of the Poisson brackets between coefficients of charac- 
teristic polynomials ^(A) (3.13) and B{X) (3.12). 

Now we have to determine an inverse transformation from variables ui,2 and pu^ 2 to initial 
variables p, M. Firstly, using the definitions (3.13), (3.12) and the relation (3.10), we express 
five linear functions (3.9) via variables of separation and Casimir functions 

- Vu2 2ti - 3ui 2ti - 3n2 

ui - U2 3(ni - U2) 3(ni - U2) 



{ui-U2)Ci ( tI\ Pm-Pu2 , 2ti Pu^U2 - PU2UI , {Puj_U2 - Pu2Ulf 
P2 = 77 T ~ "^2 + 7^7 r + -T, \- 



^Pui-Pu2) V 9y2(Mi-U2) 3 U1-U2 4:{Puj_ - Pu2){ui - U2)' 

Ul + 2U2 2UI+U2 /o-,A\ 

0-0 = Pui Pu2, (3.14) 

Ui — U2 Ui — U2 

Tipi 2{pu^Ui - Pu2U2)ti {Put, - Pu2)uiU2 

c\ = -P2 + ^7 ^ \ • 

3 3(ui - U2) u\ - U2 

It is easy to see that pk and (T^ are symmetric functions in u\^2i Pui 2- 
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Secondly, we determine the initial variables p, M as functions on and cr^,: 

Pi = Ci{ajakPo - {aj + ak)pi + P2), i = 1,2,3, {i, j, k) = (1,2,3), 

Mi = — (C2 - cri{aipo - pi) + ao(p2 - npi + ai{aj + ak)po)). (3.15) 
Po 

Here Cj are entries of the vector c (3.8), Ci^2 are the Casimir functions on e*(3) (2.3) and k) 
means the permutation of (1,2,3). 

Matrices Fi^2 (3.11) in canonical variables of separation look like 

2:>-' 

where the Stackel matrix S is equal to 

^\ 

5 = 2ti - 3ui 2ri - 3^2 . (3.I6) 



The notion of the Stackel matrix S allows us to easily get the separated relations (1.1) and prove 
that canonical variables u, p„ are the variables of separation for the Steklov-Lyapunov system. 

Proposition 6. In the Steklov-Lyapunov case the canonical variables ui^2 (3.12) andp^i 2 (3.13) 
satisfy the following separated relations 

^{u,Pu) = Hi + H2 + (f3{u)pl + (l)3{u) =0, = Ml,2, P«=P«i2) (3-17) 




,2 3 

where cubic polynomials <P3{u) and (psiu) are equal to 



^3{u) 



(U^ U 4: \ 

y + 3 (^1 - 3r2) - — (n - 3ai)(ri - 3a2)(ri - Sas) j , 

The proof consists of substituting integrals of motion Hi^2 (2-7) in terms of variables of 
separation (3.15) into the separated relations (3.17). 

So, in the Steklov-Lyapunov case equations of motion are linearized on Jacobian of the genus 
two hyperelliptic curve defined by the equation ^{u,pu) = (3.17) and the system of the 
Abel-Jacobi equations has the standard form 

du du 

+ / — ^ = /3it + 7i. 



oo p{u)^pz{u) Joo p{u)lp3{u) 
"1 udu udu „ 

i . I . + / -TT-^\ = + 72- (3.19) 

o P[U)^P2,(U) 7oo P{u)f3{u) 

Here p{u) means the function p^ on u obtained from the equation (3.17), /3i^2 are certain con- 
stants depending only on the choice of the Hamiltonian (Hi or H2) and 71^2 are two constants. 
Solving these equations with respect to 1*1,2(^)71,2) and substituting these solutions into the ex- 
pressions (3.14) and (3.15) we finally get the initial variables p, M as functions on time variable t 
and six constants -ffi 2 , C12 and 71 2 . 
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Remark 5. In order to give an explicit theta-functions solution, one can apply the standard 
machinery of the Weierstrass root functions describing inversion of the hyperelliptic quadra- 
tures (3.19), completely similar to solution of Jacobi's geodesic problem or Neumann's particular 
case of the Clebsch system [5, 8, 27]. 

Using shift of coordinates 
2 

Ui -(ai + 02 + as - 3ui) 
we can rewrite separated relations (3.17) in the following form 

2 _ Ciu^ + {Ci tr A - C2W + + 4^2 

detiu-A) ' ^"^-^"^ 

where 

= -Hi - Ci tr (A^ + 4A^) + 2C2 tr A, H2 = H2 + 2Ci det A - 2C2 tr A^ . 

Separated relations (3.20) have the same form as the Kotter separated relations (2.14). However, 
we have to point out that if p and M are real variable, then ui^2 are complex functions in contrast 
with the real Kotter variables vi^2- In terms of ui^2 and p^^ ^ symmetric functions on the Kotter 
variables f 1^2 look like 

P{ui,U2,Pui,Pu2) R{ui,U2,Pui,Pu2) 

Vl+V2 = yr-, r, V1V2 = —f V 

Q[Ui,U2,Pui,Pu2) J [Ul,U2,Pui,Pu2) 

Here P, Q, R, T are the sixth-order polynomials in momenta, Q, T are the sixth-order polyno- 
mials in coordinates, whereas P and R are the seventh- and eighth-order polynomials in coor- 
dinates. We did not find a foreseeable expressions for these polynomials or their combinations. 
In any case variables u, pu and v, p^ are related by non-trivial canonical transformation (1.3). 

Remark 6. The Steklov-Lyapunov system on e*(3) coincides with the Steklov system on so* (4) 
after some linear change of phase variables [3]. It is a twisted Poisson map, which permutes first 
and second Lie-Poisson brackets on e*(3) and so* (4) [19]. 

We suppose that variables mi,2 and Pux,2 coincide with the complex variables of separation for 
the Steklov system on so(4) introduced in [6, 13] up to this change of variables and transformation 
of the canonical momenta associated with permutation of the Poisson brackets. We thank one 
of the referees for the reference on these papers. Algebro-geometric relations of this complex 
coordinates ui^2 with the real Kotter coordinates vi^2 is discussed in [8]. 

3.3 The Rubanovsky system 

Let us consider a nontrivial integrable generalization of the Steklov-Lyapunov system discovered 
by Rubanovsky [16] 

Hi = Hi + 2{b,p), H2 = H2 + {b,{trA-A)p)-{b,M), (3.21) 

where Hi^2 are given by (2.7) and b = (61, 62; ^3) is a constant vector. This deformation describes 
the motion of a gyrostat in an ideal fluid under the action of the Archimedes torque, which arises 
when the barycenter of the gyrostat does not coincide with its volume center. The problem of 
separation of variables for the Rubanovsky systems was unsolved up until now. 

The Rubanovsky integrals of motion are non-homogeneous second-order polynomials and, 
therefore, it is natural to solve the equations (3.1) and (3.2) in the space of non- homogeneous 
second-order polynomial bivectors and in the similar space of rational bivectors. 
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Proposition 7. Integrals of motion Hi^2ip, M) (3.21) are in bi-involution (3.1) with respect to 
the Poisson brackets associated with the polynomial Poisson bivector 



Pi=P[ + {b,c) 



and the rational Poisson bivector 



/o 














\ 


















































P3 


-P2 











-P3 





Pl 


\o 








P2 


-Pl 


/ 



(3.22) 



1 



1 



{c,p) {c,p) 



, , o(l + a^(c, d)) , « 1 „ 

(c,M) + A^-,^(cxp,d))P, 



(3.23) 



where P[ is the Poisson bivector for the Steklov-Lyapunov system (3.4) in which a = 1 and 
vectors c and d are given by (3.8). 

The proof is a straightforward verification of the equations (3.1) and (3.2). 

Remark 7. If c = (ci, C2, y^— — c^) is an arbitrary vector, d = (cia~^, 0, 0) and a — t- oo, then 
the Poisson bivectors P[ 2 yield bi-Hamiltonian structures on e* (3) associated with the Lagrange 
top [20]. 

Following the same line as previously, let us calculate the control matrices associated with 
the Hamiltonians Hi^2 (3.21) and the Poisson bivectors P{ 2 (3.22), (3.23) 



(&,c) 
1 



Fo — Fo 



{b,c) fo 



2po yi 



Eigenvalues of these matrices are new canonical variables of separation ui^2 and pui 2, so that 



Fi = S 



2p 



2p2, 



Fo = S 



ui 

U2 



s- 



where S is the same Stackel matrix (3.16). These variables are simply related with the previous 
one 



Ul,2 



Pui : 




{Ui - U2)' 



4(&,c) 



Initial variables p, M are the same functions on pk and ak 

+ afc)/Oi +/02), « = 1,2,3, (i, j,/c) = (1,2,3), 



Pi = Ci[ajakPo - 

-^{C2 - cri{aiPQ - Pl) + cfq[p2 - npi + ai{aj + ak)po)), 
Po 



where four functions Pk{u,Pu) = Pkiu,p) and ao{u,pu) = ao{u,pu) is given by (3.14): 



Po 

P2 



q P"! Pu2 

Ul — U2 

Ui + 2U2 ^ 

—Pui 

Ul - U2 

(ui -U2)Ci 



Pl 



2ti — 3ui 



3{ui - ^2)^"^ 3(ni - ^2) 



2ti — 3u2 



P 



U2 ) 



2ui + U2 . 

~ ; Pu2 5 

Ul - U2 



4(p 



Ul 



p 



■U2) 



T2 + 



Pui - Pu2 ^ '^TiPuiU2 -PU2U1 _^ 



9 / 2 Ml - U2 



Ul — U2 



{PuiU2 -Pu2Ulf 
4(Pwi -Pu2){ui - U2) 
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and one function is different 

{b,c) . Tipi 2{pu^Ui - Pu2U2)ti {Pui - Pu2)uiU2 

c^l = P2 + — :7Ti \ : : ■ 

Z 3 3(^1 — U2) Ul — U2 

This shift of ai acts only on M variables (3.15). 

Proposition 8. In the Rubanovski case the separated relations have the following form 

Hu,Pu) =i^--j-jHi+H2 + Mu)pI - + ^-^^ + Q J Pu (3.24) 

, /-N (b.c) f u Ti\ (b,c)(b,Ac) . . . . 

+ 03(n)-^ (^- + ^J +^^^^^ ^=0, ^X = ni,2, Pu=Pu,,2- 

Here cubic polynomials 993 and (p^ are given by (3.18) and D = tr A — 3 A. 

The proof consists of substituting integrals of motion Hi^2 in terms of variables of separa- 
tion u, Pu into the separated relations (3.24). 

As above, we can prove that the equations of motion are linearized on the Jacobian variety 
of the genus two hyperelliptic curve defined by (3.24). For the brevity, here we omit the explicit 
expressions for a base of holomorphic differentials and the corresponding Abel-Jacobi equations, 
which may be easily obtained using modern computer algebra software. 

Calculation of the separating variables and the corresponding algebraic curve for the Ru- 
banovski gyrostat is a new result, which allows us to make conclusion an applicability of bi- 
Hamiltonian methods to study finite-dimensional integrable Hamiltonian systems. 
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